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'^ . Abstract 

'■^ \ We investigate the role of nuclear spin in elastic scattering of Dark 
'^ ■ Matter (DM) neutralinos from nuclei in the framework of the Mini- 
ng . mal SUSY standard model (MSSM). The relative contribution of spin- 

^ \ dependent axial-vector and spin-independent scalar interactions to the 

event rate in a DM detector has been analyzed for various nuclei. Within 

^vg . general assumptions about the nuclear and nucleon structure we find 

^— I \ that for nuclei with atomic weights A > 50 the spin-independent part 

^-r ' of the event rate Rsi is larger than the spin-dependent one Rg^ in the 

Q"^ ■ domain of the MSSM parameter space allowed by the known experi- 

"^ \ mental data and where the additional constraint for the total event rate 

Q^' R = Rsd + Rsi > 0.01 is satisfied. The latter reflects realistic sensi- 






Dh! 



^ 



X 



tivities of present and near future DM detectors. Therefore we expect 



D , equal chances for discovering the DM event either with spin- zero or with 



spin-non-zero isotopes if their atomic weights are Ai ~ ^2 > 50. 



^ This work was supported in part by the Russian Foundation for Fundamental Research 
(93-02-3744). 
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1 Introduction 

Analysis of the data on distribution and motion of astronomical objects within 
our galaxy and far beyond indicates presence of a large amount of non- luminous 
dark matter (DM). According to estimations, it constitutes more than 90% 
of the total mass of the universe if a mass density p of the universe close 
to the critical value pcrit is assumed. The exact equality Vt = p/pcrit = 1, 
corresponding to a flat universe, is supported by naturalness arguments and 
by inflation scenarios. Also, in our galaxy most of the mass should be in a 
dark halo. Detailed models predict a spherical form for the galaxy halo and a 
Maxwellian distribution for DM particle velocities in the galactic frame. The 
mass density of DM in the Solar system should be about p ~ 0.3GeV-cm~^ 
and the DM particles should arrive at the earth's surface with mean velocities 
V ~ 320 km/sec, producing a substantial flux $ = p ■ v/M ($ > lO'^cm"^ 
sec~^ for the particle mass M~ 1 GeV). Therefore one may hope to detect DM 
particles directly, for instance through the elastic scattering from nuclei inside 
a detector. 

The theory of primordial nucleosynthesis restricts the amount of baryonic 
matter in the universe to 10%. Thus a dominant component of DM is non- 
baryonic. The recent data by the COBE satellite [T[] on anisotropy in the 
cosmic background radiation and the theory of the formation of large scale 
structures of the universe lead to the conclusion that non-baryonic DM itself 
consists of a dominant (70%) "cold" DM (CDM) and smaller (30%) "hot" DM 
(HDM) component H, iJ. 

The neutralino (x) is a favorable candidate for CDM. This is a Majorana 
ix'^ = x) spin-half particle predicted by supersymmetric (SUSY) models. 

There are four neutralinos in the minimal SUSY extension of the standard 
model (MSSM). They are a mixture of gauginos (W3, B) and Higgsinos (-^1,2), 
which are SUSY partners of gauge {W^, B) and Higgs (-^1,2) bosons. The 
DM neutralino x is the lightest of them. Moreover, x is assumed to be the 
lightest SUSY particle (LSP) which is stable in SUSY models with i?-parity 
conservation. 

The problem of direct detection of the DM neutralino x via elastic scat- 
tering off nuclei has been considered by many authors and remains a field of 



great experimental and theoretical activity [H-[|TT|. 

The final goal of theoretical calculations in this problem is the event rate R 
for elastic x-nucleus scattering. In general, the spin-dependent {Rsd) and spin- 
independent {Rsi) neutralino-nucleus interactions contribute to the event rate: 



R = Rsd + Rsi- Rsd vanishes for spinless nuclei and this fact is often regarded 
as a reason to assert spinless nuclei to be irrelevant for the DM neutralino 
detection as giving a much smaller event rate. One can meet this statement 
in the literature. However, this is right only if the spin-dependent interaction 
dominates in elastic neutralino scattering off nuclei with non-zero spin. 

In this paper we address the question on the role of nuclear spin in the 
DM neutralino detection. We investigate this problem in the framework of 
the MSSM. We avoid using specific nuclear and nucleon structure models but 
rather base our consideration on the known experimental data about nuclei 
and nucleon. It allows us to free the consideration of theoretical uncertainties 
specific for the structure models. To restrict the MSSM parameter space we 
use experimental constraints on SUSY-particle masses, the cosmological bound 
on neutralino relic abundance and the proton life-time constraint. 

We have found that Rgi contribution dominates in the total event rate R 
for nuclei with atomic weight A > 50 in the region of the MSSM parameter 
space where R = Rsd + Rsi < 0.01. The lower bound 0.01 is far below the 
sensitivity of realistic present and near future DM detectors. Therefore we can 
exclude the region where R < 0.01 as invisible for these detectors. 

We do not expect a crucial dependence of the DM event rate on the 
nuclear spin for detectors with target nuclei having an atomic weight larger 
than 50. As a result, we expect equal chances for J = and J 7^ detectors 
to discover DM events. In particular, this conclusion supports the idea that 
presently operating /5/3-detectors with spinless nuclear target material can be 
successfully used for DM neutralino search. These highly developed set-ups 
(for a review see fl^), operating under extremely low background conditions, 
use detection technology which is suitable for the DM search. 

2 General Properties of the Neutralino - Nu- 
cleus Interactions. 

A DM event is elastic neutralino-nucleus scattering causes the nuclear recoil 
detected by a detector. The event rate per unit mass of the target material 
depends on the distribution of the DM neutralinos in the solar vicinity and 
the cross section o"e/(x^) of neutralino-nucleus elastic scattering. One can 
calculate aei{x^) starting from the neutralino-quark effective Lagrangian. In 



the most general form it can be given by the formula 

Leff = ^{^(g) - A{q)} ■ X7/.75X ■ 97^75? + ^T" ' ^i ' XX ■ QQ, (1) 

q J-Vlw 

where terms with the vector and pseudoscalar quark currents are omitted be- 
ing neghgible in the case of the non-relativistic DM neutrahno with typical 
velocities v^ ~ 10~^c. Following the tradition we retain in the first term the 
difference of two coefficients representing just one independent parameter of 
the Lagrangian. The coefficients B{q),A{q),C{q) depend on the SUSY model 
and will be considered in the next section. Here we survey general properties 
of X"^ scattering following from the Lagrangian (|I]). 

To calculate a^iix^) o^^ should average the x-l interactions sequentially 
over the nucleon and the nuclear structure. The first and the second terms 
in Leff ([5) averaged over the nucleon states give the spin-dependent and the 
spin-independent matrix elements M^^^ and M^j, respectively. For the spin- 
dependent amplitude we have 0, |]§|: 

m,, = ASj,^n) E {B{q)-A{q)}Aq, (2) 

<jep(n) 

where S-^ and Sp(n) are the neutralino and proton (neutron) spin operators; 
Ag are the fractions of the nucleon spin carried by the quark q. The standard 
definition is 

< P{n)\q^>'-i,q\p{n) >= 2S^^^^Aq, (3) 

where Sp,^^ = (0,5'p(„)) is the 4-spin of the nucleon. The parameters Ag (for 



the proton) can be extracted from the EMC ||T3[ and hyperon data [|T4| : 



Am = 0.77 ±0.08, Ad = -0A9 ±0.08, As = -0.15 ± 0.08. (4) 

The relevant values for the neutron can be found from (^) by the isospin 
symmetry substitution u ^ d, d ^ u. 

The spin- independent matrix element has the form 0], [QQ 



M 



SI 



ruu + md 



+ ^{l - f - f){C{c) + C{h) + C{t)) 



MpM. . yy . ^^ 



w 



"'^When this paper had been completed we received a paper ref. llQ] with more refined 
treatment of the spin- independent matrix element. 



where the parameters / and / are defined as follows: 

<p{n)\{mu + md)iuu + dd)\p{n)> = 2/Mp(„)^^, (6) 

< p{n)\msss\p{n) > = /Mp(„)^^. 



The values extracted from the data are |T5|,[|I^: / = 0.05 and / = 0.2. 

Averaging @, (|^) over the nuclear states |A > we deal with the following 
matrix elements at vanishing momentum transfer: 

< A|Mp(„)M|A > = MaAA, (7) 

< A\Spi^ri)\A > = X<A\J\A>. 



Here J is the nuclear spin. On the basis of the odd-group shell model |T^ 
(essentially somewhat relaxed single particle shell model) the parameter A can 
be related to the nuclear magnetic moment, fi, as follows 

g' - 9^ 

where g'' = 1(0) and g'^ = 5.586(— 3.826) are orbital and spin proton (neutron) 
^f-factors. Then one can extract values of A for various nuclei from the experi- 
mental data on nuclear magnetic moments f\ We use in this paper the values 
of A as presented in ref . [0] . 

For large M^ and Ma the momentum transfer may be comparable to the 
inverse radius of a nucleus and then we have to take into account the finite size 
effect (see also [0). It can be done by introducing the coherence loss factor 

m- 

0.573/ exp(-^)erf(y5^) 
where 






9 {M^ + MaY' 

Here cr^ is the dispersion of the Maxwellian neutralino velocity distribution a = 
0.9 ■ 10~^. To obtain the coherence loss factor for spin-independent scattering 



^ A more direct way of calculation based on the theory of finite Fermi systems is presented 
inrtisl. 



we take r 
AM: 



r charge in (|), wheie r charge IS the ruis charge radius of the nucleus 



Tchargc = (0.3 + 0.89M'f) fill 



(10) 



The coherence loss factor for spin-dependent scattering is given by 



with r 



spin- 



The rms spin radius of the nucleus A can be estimated as 



^spin = ^ ■ r charge with ^ ?^ 1.25 froui harmonic well potential calculations P|. 

Finally we arrive at the formula for the event rate of elastic neutralino- 
nucleus scattering in the detector per day per unit mass of the target material: 



R 



J^si ~r J^sdi 



where the spin-dependent and spin-independent parts are: 
Rsd = 5.8-10i°-AV(J+l)C(r,pi„);W 



i?o 



1.44 ■ 10 
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AV(J+l)C(r,p,„ 



2 .T) 
sd ^y 



^M, 



w 



Circharge) " M^ ■ V. 



(11) 
(12) 

(13) 



The common kinematic factor V and properly normalized nucleon matrix 
elements J^si^M.sd are defined as: 



V 



AM^Ma 



< \ve\ > 



events 



M 



(M^ + M^)2J L.3G'e1/ ■ cm-3J V^2Qkm/s^ kg ■ day 



XX ■ ^^, 



Mw 

A- Msd- SxSp(n) 



(14) 

(15) 
(16) 



For the definition of M,,-, M.^ see formulae 



Here p ~ 0.3 GeV-cni 



-3 



is the DM neutralino density in the solar vicinity and < \ve\ >~ 320 km/s is 
DM neutralino averaged velocity at the earth's surface. 

To study the role of nuclear spin in elastic x-nucleus scattering we introduce 
the ratio 

K = Rsd/Rsi (17) 

characterizing the relative contribution of spin-dependent and spin-independent 
interactions. From the practical point of view it determines the expected rel- 
ative sensitivity of DM detectors with spin-non-zero (J 7^ 0) and spin-zero (J 
= 0) nuclei as target material. li k < 1, then detectors with spin- non- zero and 
spin-zero target materials have approximately equal sensitivities to the DM 
signal, whereas if k > 1 then, the spin-non-zero detectors are more sensitive 
than the spin-zero ones. 



Let us consider separately the dependence of n on the nuclear structure 
and the choice of a specific SUSY model. We may write: 



^ VAVsusy^ 



where 



C{rspin)M\ 



w 



Here rjA is a factor depending on the properties of the nucleus A; r]'f^sy is defined 
by the SUSY-model which specifies the neutralino composition and the inter- 
actions with matter. The SUSY-factor also depends on the shell-model class 
to which nucleus A belongs, being tj'^usy ^^ ^^e shell-model "neutron" (^He, 
29Si, ^^Ge,...) and ifsusv for the shell-model "proton" (^^F, ^sq, ^osti,...). 

Fig.l depicts the nuclear factor tja versus the atomic weight A. The error 
bars represent the interval of the neutralino masses 20 GeV< M^ < 200 GeV. 



The lower bound corresponds to the present experimental constraints p0|-p2 



The upper bound is taken to include recent estimations for the mass of the 



cosmologically favorable neutralino [23|. It follows from fig.l that ?7^ < 1 for 



A > 50. Thus at A > 50 there is no nuclear structure enhancement of the 
spin-dependent event rate as compared to the spin-independent one. 
The next is an estimation of the SUSY-factor r^f^"^. 

3 Specific SUSY-model Predictions. 



To estimate the factor tjsusy in (|20| ) one should calculate the parameters A{q), 
B{q) and C{q) of the effective Lagrangian (|l]) in the specific SUSY model. We 
will follow the MSSM. This model is specified by the superpotential and "soft" 



SUSY breaking terms |24 



The effective low energy superpotential is: 

W= J2 {huH2QU + hDH^QD + hLH^LE)+fiH^H2. (21) 

generations 

Hi and H2 are the Higgs fields with a weak hypercharge Y = —1, +1 respec- 
tively. 



SUSY breaking in the "hidden" sector of N=l supergravity produces "soft'' 
supersymmetry breaking terms in the scalar potentiah 



ysoft = E ^^i\^i\^ + huAuH^QU + hDAoH^Qb + hLALH^LE 

i=scalars 



:22i 



+I1BH1H2 + h.c. 



and a "soft" gaugino mass term 

1 
"2 



C 



FM 



M^BB + M:^W^W^ + Mg^"^" 



— h.c. 



(23) 



The model is also characterized by the set of boundary conditions at the 
unification scale Mx'- 

Ajj = Ad = Al = Ao, (24) 

(25) 



m 



Hi 



m 



H2 



m, 



m, 



Ulr 



m, 



m^ = mo, 



M^ = M2 = M3 



m 



1/2? 



(26) 
g,{Mx) = g^iMx) = 92{Mx) = Qgut, (27) 

where 9^,92,91 are the SU{3) x SU{2) x t/(l) gauge coupling constants equal 
to Qgut ^^ ^^6 unification scale Mx- At the Fermi scale Q ~ M^ these 
parameters can be evaluated on the basis of the 1-loop renormalization group 
equations (RGB) ^,^. 



The neutralino mass matrix in this model has the form |p4i: 



/ 



M. 



M2 


-MzCyirSp 

V MzCwCp 





-MzSwCp 



-MzCy^sp Mzc^cp \ 

MzS^Sp -MzC^rSp 

-fi 

-/i y 



:28) 



where c^r = cos^vv, s^y = sin 9w, t^r = tan6'iy, sp = sin/3, c/j = cos/?. The 
matrix is written in the basis of fields (W^, B, H2, H^). As usual, M2, fi are 
the gaugino mass and the Higgs mixing parameter; the angle f3 is defined by 
the vacuum expectation values of the neutral components of the Higgs fields: 
tan/? =< H2 > / < H^ >= f2/'yi. Diagonalizing the mass matrix ( pH]) we 
obtain the lightest neutralino of the mass M^ with the field content 



X = aW^ + f3B + 7^2° 



6Hl 



To calculate the low-energy neutralino-quark interactions we also need to have 
the spectrum of squarks q and Higgs particles at the Fermi scale. Their mass 
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matrices depend on the soft SUSY breaking parameters. We obtain them from 
the boundary conditions at the GUT scale Mx (PD- (|26|) as a solution of the 
1-loop RGE. 



We analyze the Higgs sector of the MSSM at the 1-loop level [^. In 
the analysis we take into account i^ — ifi, hi — Br mixing between the third- 
generation squarks. Diagonalization of the Higgs mass matrix gives three neu- 
tral mass-eigenstates. There are two CP-even states, H, h, with the masses 
''TT'H, fTT'h and the relevant mixing angle an and one CP-odd A with the mass 
rriA- We take the mass rriA as an independent free parameter of the MSSM. 

A complete list of essential free parameters of the MSSM is 

t&iaP,Ao,B,fi,mi/2,mA,mo,mt. (29) 

Having a particle spectrum one can derive the effective Lagrangian L^ff of low- 
energy neutralino-quark interactions. As discussed in the previous section, its 
general form is given by eqn. (|1|). In the MSSM the first term of L^ff is 



induced by the Z-boson and q exchange ||28| whereas the second one is due to 



the Higgs particle [^ and q exchange as well as qi — qn mixing ||^,pO 
The coefficients of L^ff are 

|7|2_|5|2 Yl 

A{q) = TTT2 — {gi sm 9w + 92 COS 9w) i^Y 3-^ ^^^^w~ '^^92 COS 9 w) 

l H2T, + (3g,^\^ 1 ml M Ts j - Ts 

+ 2 ml,-M-^ +2m?^-M2^ v^ '^'+ vf '^' J' ^^"^ 

|7|2-W|2 Yr 

B{q) = ^^ — (5fism6'iy + 5f2Cos6'i^)(— 5fism6'H/) (31) 

1 M^gdf? 1 m\ i + Ts 1-^3 , |2. 

2m?^-M2 2m?^-M2^ vl 1^' + vl I'' J' 



^2 7T.r/l , r^-^cosan ,1 sina^^ 



C{,) = T^F,[it + T,r-^^-it-T,r-!^] (32) 



Here 



g2 ag2T3 + ^Pgi ^giP 1 ^.7 .1 . ^ . 

4^ m?^-M2 m?^-M2J^4 '^sin/? 4 '^'cos/3^ 



Fq = {a — P ta.n 9^) (7 cos an + S sin an)- (33) 
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In these formulae we ignore ql — qr mixing since they give a small contribution 
according to the estimation of ref . . 

The procedure we use for the neutralino mass matrix diagonalization always 
leads to positive mass eigenvalues and to either real or pure imaginary values 
of the coefficients a, P, 7, 6. Therefore in formulae (|30D, (|3l|) the absolute values 
of these coefficients appear. 

Now we are ready to calculate the rjsusy-i^ctoi (|3) substituting the defini- 
tions (ig)-(|3|) in formula (|g). 

To get complete information about possible values of the T^s^sy-factor we 
scan the MSSM parameter space within the constraints imposed by the exper- 
imental data and some general theoretical principles. The well known experi- 
mental constraints 



are summarized in the Table. 



Table: Present Limits on Supersymmetric Particles. 

(Table is taken from ref. ||31||.) 





Bound on 




Particle 


Particle Mass 
(GeV) 


Source 


Xl 


18.4 


f Based on the LEP non-observation of x i 


xl 


45 


I and Xl, CDF non-observation of g, and 


Xl 


70 


*- the assumption of gaugino mass unifica- 


v° 


108 


tion. 


45.2 


LEP 


xt 


99 


See neutralino mass limits above. 


V 


41 


Assumes the z/ decays are invisible (other- 
wise M i> < 32 GeV). Based on LEP mea- 
surement of r{Z —>■ invisible final states). 


e 


45 


LEP; assumes M^o < 41 GeV 


f^ 


45 


LEP; assumes M^o < 41 GeV 


f 


45 


LEP; assumes M^o < 38 GeV 


<? 


45 


LEP; assumes M^o < 20 GeV 




74 


UA2(anyMg) 




-95 


CDF (M,- < Mg) 


9 


79 


UA2 {M-g < Mg) 




-95 


CDF 



The constraints are given for masses of squarks q, the gluino g, charginos 
X"^, the neutralino x, charged sleptons /^, the sneutrino z/, lightest CP-even h 
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and CP-odd A Higgs bosons. We also include the constraints 

1.12 < tan/? < 4.7. (34) 

The lower limit follows from the finiteness condition for the top Yukawa cou- 
pling Yt. If 

sml3> {mt/200GeV), (35) 

then Yt is finite up to the unification scale Mx- For rrit = 150 GeV we get the 
lower limit in (B5). The upper limit in this formula is expected from proton 



stability considerations |32]. 



From the "naturalness" arguments [^ we may choose 



mq,mg < 1 TeV, (36) 

where rrif is the mass of any sfermion /. The choice of the interval for the 
neutralino mass 

20 GeV <M^< 200 GeV 

was already explained at the end of section 2. 

The additional constraint we use in the analysis of the role of nuclear spin 
is the constraint on the realistic sensitivity of the DM detector. In terms of 
the total event rate R we choose the sensitivity to be not better than: 

R > 0.01 /"^f . (37) 

kg ■ days 

We do not expect the DM detectors to go below this lower bound in near 
future. Therefore the constraint (|37|) reflects the realistic capacities of the 
present and near-future set-ups. It excludes the region in the parameter space 
corresponding to the low-level DM signals inaccessible to these detectors. 

We have performed a numerical analysis of the MSSM parameter space 
within the above-defined constraints. In fig. 2 the typical behavior of the rjsusy- 
factor in particular domains of the MSSM parameter space is presented. The 
following upper bound for the SUSY-factor in eqn. (llHD was found: 



risusy < 1.2. (38) 

Combining this result with the values of the nuclear factor r^^ represented in 
fig.l we conclude that 



K = Rsd/Rsi = VAVsit < 1 for nuclei with A > 50 (39) 
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at a detector sensitivity up to R > 0.01. The tendency is that at higher 
sensitivities (lower R accessible) we get k < 1 for heavier nuclei. 

As a by-product of our analysis in fig. 3 we also give the event rate for some 
nuclei of special interest in DM search. 

We do not take into account possible rescaling of the local neutralino den- 
sity p which may occur in the region of the MSSM parameter space where 
Qh"^ < 0.05 0. This effect, if it took place, would essentially reduce the event 



rate R [TT|. Of course, it has no influence on the ratio k in the formula (|T] 
and on our conclusion about the role of nuclear spin. Plots in fig. 3 correspond 
to a situation when neutralinos constitute a dominant component of the DM 
halo of our galaxy with the density p = 0.3 GeV-cm~^ in the solar vicinity. 

4 Conclusions 

In the framework of general assumptions about the nuclear and nucleon struc- 
ture considering the MSSM as the basis for description of the neutralino prop- 
erties we have drawn the following basic conclusions. 

For sufficiently heavy nuclei with atomic weights A > 50 the spin-independ- 
ent event rate Rsi is larger than the spin-dependent one Rsd if low-level signals 
with total event rates R = Rsd + Rsi < 0.01 are ignored. This cut condi- 
tion reflects the realistic sensitivities of the present and the near-future DM 
detectors. 

The main practical issue is that two different DM detectors with (J = 0, Ai) 
and with (J 7^ 0, A2) nuclei as target material have equal chances to discover 
the DM event if Ai ~ A2 > 50. 

Another aspect of the DM search is the investigation of the SUSY-model 
parameter space from nonobservation of DM events. Apparently, in this case 
experiments with J 7^ nuclei are important since they provide new informa- 
tion about the SUSY model parameters from Rg^ which is inaccessible in J = 
experiments. 

The results presented above were obtained in a specific SUSY-model. There- 
fore it is a natural question whether our basic conclusions hold for other popu- 
lar SUSY-models. We plan to investigate this question in a subsequent paper. 
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Figure Caption 



Fig.l The nuclear factor rjA versus the atomic weight A. The error bars repre- 
sent the interval of the neutralino masses 20 GeV< M^ < 200 GeV. 

Fig. 2 The SUSY factor rjsusy versus neutralino mass M^ at various values of the 
MSSM free parameters, (a) and (b) for nuclei with proton and neutron 
shell model structure, respectively. 

Fig. 3 Samples of plots for the total event rate for neutralino elastic scattering 
off (a) equal parts mixture of "^^Ge +^^ Ge and (b) ^^F nuclei. 
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